Reducing Polarization Mode Dispersion With Controlled Polarization Rotations 
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One of the fundamental limitations to high bit rate, long distance, telecommunication in optical 
fibers is Polarization Mode Dispersion (PMD). Here we introduce a conceptually new method to 
reduce PMD in optical fibers by carrying out controlled rotations of polarization at predetermined 
locations along the fiber. The distance between these controlled polarization rotations must be less 
than both the beat length and the mode coupling length of the fiber. This method can also be 
combined with the method in which the fiber is spun while it drawn. The incidence of imperfections 
on the efficiency of the method is analysed. 

PACS numbers: 42.81. Gs 42.81. Bm 



Introduction. Polarization Mode Dispersion (PMD) 
arises when two orthogonal states of polarization prop- 
agate at different velocities in optical fibers. PMD is 
one of the fundamental limitations on high speed, high 
bit rate communication in fiber systems because it dis- 
torts the shape of light pulses, and in particular induces 
pulse spreading. Hence much effort has gone into reduc- 
ing PMD in optical fibers, see 0,0 for reviews. 

PMD arises because of uncontrolled stresses or 
anisotropies induced in the fiber during the manufac- 
turing process and during deployement. These cause 
unwanted birefringence, and hence PMD. This residual 
birefringence changes randomly along the fiber, resulting 
in random mode coupling as the light propagates along 
the fiber. Because of this statistical process, the effects 
of PMD -such as pulse spreading- increase as the square 
root of the propagation distance]^, 0, S IS ■ This sta- 
tistical process makes it extremely difficult to correct the 
effects of PMD after the light has propagated through a 
long length of fiber. For this reason one rather tries to 
reduce the PMD of the fiber itself. 

Two main methods have been devised to reduce PMD 
in optical fibers. The first is to minimize asymmetries 
in the index profile and stress profile of the fiber. To 
this end the manufacturing process has been steadily im- 
proved. The second method is to spin the fiber during 
the manufacturing process as described in 0, . Spin- 
ning does not reduce the anisotropies in the fiber. Rather 
it modifies the orientation of the anisotropies along the 
fiber in such a way that after one spin period the effec- 
tive birefringence and PMD is reduced. The effects of 
spinning on PMD were analyzed in e.g. [ToL fill fl2L IT^|. 
Residual birefringence can be characterized by the fiber's 
beat length, which in present day fibers is of order 10m. 

In the present work we introduce a conceptually differ- 
ent method for reducing PMD in optical fibers. Our basic 
idea is to introduce controlled polarization rotations at 
predetermined locations along the fiber in such a way 
that the effects of PMD are reduced. This is paradox- 



ical since the aim is to reduce birefringence and PMD, 
and we claim to achieve this by introducing additional 
polarization rotations. Our main idea is that after sev- 
eral such controlled polarisation rotations, the state of 
polarisation of the light has been oriented along many 
different directions on the Bloch sphere in such a way 
that the effects of the fiber birefringence averages out. 
We will show that this can be a highly effective method. 
It can also be combined with the method of spinning the 
fiber. 

We will not be concerned here with how the controlled 
polarization rotations are implemented in practice. We 
note that there are many ways of inducing controlled 
birefringence in fibers, and hence controlled polarization 
rotations. Typically these imply inducing stress in the 
fiber, for instance through compression, twisting, etc. We 
leave for future work to determine which implementation 
of controlled polarization rotations is most suitable for 
reducing PMD. 

The method presented here for reducing PMD in op- 
tical fibers is inspired by the pulse sequences (e.g. spin 
echo) that have been developed in Nuclear Magnetic Res- 
onance (NMR) to reduce the effects of imperfections, see 
for instance It is also related to the "bang-bang" 

technique introduced by Viola and Lloyd to reduce de- 
coherence in the context of quantum information |lfil flfij | . 
Wu and Lidar [J^ suggested applying bang-bang control 
to reduce losses in optical fibers through controlled phase 
shifts, but their proposal is impractical because the con- 
trol operations need to be implemented over distances 
comparable to the wavelength of light. On the other 
hand in the present case the control operations need to 
be implemented over distances of order lm, which makes 
reducing PMD with controlled polarisation rotations a 
practical proposal. An important conceptual difference 
which we stress is that in the bang-bang technique one 
wishes to average to zero the coupling between the sys- 
tem and the environment, whereas in the present case it 
is the coupling between the polarisation and frequency 
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of light itself, ie. two degrees of freedom of the same 
system, which one wishes to average to zero. 

Polarized light propagating in optical fibers. 

Consider a light pulse propagating along a birefringent 
optical fiber. The pulse is centered on frequency f2 and 
has wave number K . The distance along the fiber is de- 
noted I. Its amplitude can be written as 

A(l,t)e- int+lKl (1) 

where the slowly varying envelope of the pulse A(l,t) — 

( ) is a two component vector, called the Jones 

V A 2{i,t) J 

vector, which describes the polarisation state of the light. 
We introduce the variable t' = t — l/v g where v g is the 
group velocity of the pulse, whereupon A obeys the evo- 
lution equation 

id l A = B (l)A + iB 1 (l)d tl A (2) 

where -Bo and B± are traceless, hermitian matrices de- 
scribing the birefringence and PMD along the fiber. (The 
parts of So an d B\ proportional to the identity are in- 
corporated into the wave vector K and the average group 
velocity v g . For simplicity we neglect loss, dispersion, and 
non linearities.) The term proportional to Bq describes 
the phase birefringence. Its magnitude is generally mea- 
sured by the beat length is = n/\Bo\ which is the length 
after which the polarization is brought back to its origi- 
nal state. If only Bo was present in the fiber, the state of 
polarization would evolve in an unknown way, but there 
would be no pulse spreading. The PMD is encoded in the 
group birefringence B\ which describes the difference in 
group velocities of the two orthogonal states of polariza- 
tion. This term will give rise to spreading of pulses. The 
birefringence changes along the fiber, hence both Bo(l) 
and B\{1) depend on the position I along the fiber. In 
present day low birefringence fibers both L b and the dis- 
tance over which Bq and B\ change are comparable and 
of order 10m. 

We can rewrite these matrices as 

B (l) = B (l) ■ a and B^l) = Bx(l) ■ a (3) 

where 0-,= (J J),^= (° ~q), = (o -1 ) 
are the Pauli matrices, and x, y, z denote orthogonal di- 
rections on the Poincar sphere. Denoting by A(l,u) the 
Fourier transform of A(l,t') with respect to t' , eq. 
becomes 

idiA(l,oj) = B(1,cj) ■ 5A(l,u) , (4) 
where B(l, u) = B (l) + uB^l) (5) 

is the sum of the birefringence and PMD in the fiber. 
The solution of eq. (@J can be expressed as 



with U(l) a 2 x 2 unitary matrix, called the Jones matrix 
which describes the evolution of the state of polarization 
after propagating a distance I 

Basic Principle of the method. For simplicity we 
will first suppose in the following paragraphs that B(lo) 
is independent of I, whereupon the Jones matrix is 

U(l) = e-^M-* (7) 

Let us now show that, even though we do not know 
B(1,lu), it is possible to compensate for the evolution 
so that after compensation the polarization comes back 
to its initial state. The following procedure consisting of 
four basic steps that are then repeated, realises this: 

1. The light propagates over distance I. The distance 
I is taken short enough such that \B\l << 1, hence 
the evolution can be well approximated by the first- 
order approximation: 

U(l) pa 1 - iBal = 1 - il(B x cr x + ByO-y + B z a z ) . 

2. We interrupt the evolution by flipping the polariza- 
tion around the x axis. We then let the light evolve 
over a new distance I and finally we flip again the 
polarization around the x axis. The evolution in 
step 2 is thus described by a x ll(l)a x . To first order 
in \B\l, one finds 

o x U(l)a x ps a x (l-iBal)a x = l-il{B x a x - B y a y — B z a z ) 

which effectively compensates the evolution due to 
the By and B z components. The evolution due to 
the B x component is not yet compensated; this will 
be accomplished during the next two steps. 

3. The same as Step 2, but the spin is flipped around 
the y axis. Step 3 is thus described by a y U(l)a y . 
By expanding U (I) to first order we obtain 

a y U(l)a y ps a y (l—iBal)a y = \—il{—B x a x +B y (jy—B z a z ) 

4. The same as Step 2, but the spin is flipped around 
the z axis. Step 4 is thus described by cr z U(l)a z . 
To first order we obtain 

a z U(l)a z ps a z (l-iBal)a z = l-il(-B x a x -B y a y +B z a z ) 

Putting all together, the time evolution over the four 
steps is 

U(l seq ) = a z U(l)a z a y U(l)a y a x U(l)a x U(l) (8) 
= l + 0(B 2 f seq ) (9) 

where we denote by l seq = 4/ the length of the sequence. 
The evolution is therefore effectively stopped. The proce- 
dure is then repeated again and again. Note that because 
Cy&z = i&x, eq. JBJ simplifies to 



A(l,u) = U{l i u)A(Q,w) (6) 



U(l seq ) - a z U(l)<j x U(l)a z U(l)a x U(l) . (10) 
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The method relies on the fact that the interaction, al- 
though unknown to us, does not change along the fiber 
(at least for the short lengths l seq we are considering 
here). The unknown birefringence B that is responsi- 
ble for the rotation of the polarization in the first place 
is there all the time, and affects the polarization after the 
x, y and z rotations, and brings it back to its initial state 
at the end of the sequence. 

In optical fibers the phase birefringence is of the same 
order as the group birefringence. This implies that 
Bq/VL ~ Bi. Since u>, the frequency spread of the 
pulse, is much smaller than f2, the carrier frequency, 
we can expand eq. 10 to first order in ljB\ to obtain 
U(l seq ) = 1 + 0{Bll 2 seq ) + O(uB x B l 2 seq ), Consider now 
the evolution over a fixed length L > l seq . It is given by 
U(L) = U(l seq ) L / l >°« ~ l + O(B 2 ? seg i) + O( W B 1 B / se? L) 
whereas in the absence of the control sequence it would 
have been given by 1 + BqL + ujB\L. Hence there is a 
reduction of PMD by a factor B l seq = Trl seq /LB- 

Note that as shown in it is possible to devise se- 
quences that cancel not only the terms of order Bl, but 
higher order terms as well. For instance the sequence 

U 2 (l seq ) = U(l)a x U(l)a z U(l)a x U(l)U(l)a x U(l)a z x 

xU(l)a x U(l) = l + 0(BHl q ) (11) 

cancels the evolution up to order B 3 l 3 . 

Combining Controlled Polarization Rotations 
and Spun Fibers. It is possible to reduce PMD in 
optical fibers by combining the methods of controlled po- 
larization rotations and of spinning the fiber while it is 
drawn. Here we first show how to include both methods 
in a unified description, and then we discuss how they 
can be used in combination. 

For spun fibers the evolution equation is no longer 
given by eq. J2J, but by 



id l A=Ul pin (l)B(u;)U spin (l)A 



where 



TT ■ (1\ — p^ ia ^ a y 

^ spin \ l } — c 



(12) 



(13) 



describes the spin of the fiber, a(l) is the spin function 
(the angle by which the fiber is spun), and a y is the Pauli 
matrix that generate polarization rotations around the 
circular polarization axis. To first order in B the evolu- 
tion is given by U(l) = 1 +i^'TJl# n {l')B(u))U lvin (l') 1 



from which one deduces, see |12j . that in order to reduce 
PMD the spin function must obey 

pL pi pL pi 

I dl co$[ 2a(l')dl'] = and / d/sin[/ 2a{l')dl'\ = 
Jo Jo Jo Jo 

where L is the spin period (the period of a(l)). Note that 
spinning the fiber alone can reduce linear birefringence, 
but not circular birefringence (proportional to a y ). 



If we incorporate controlled polarization rotations, we 
get the evolution equation: 

id t A = (u! pm (l)B(tu)U spm (l) + B C {1)) A (14) 

where B c (l) (proportional to Dirac delta functions) de- 
scribes the controlled polarization rotations. It is conve- 
nient to rewrite this equation in the frame that rotates 
with the spin. To this end we define A(l) = U sp i n (l)A(l). 
The rotated Jones vector obeys the evolution equation 

idiA = (b{uo) + dia(l)a y + U spin (l)B c (lpl pin (lfj A . 

^ (15) 

Thus in the rotating frame, spin is formally identical to 
a continuous controlled polarization rotation given by 
dia(l)a y , whereas the controlled polarization rotations 
B c (l) are rotated and take the form U S pi n {l)B c (l)U] pin {l). 
Therefore even though they have completely different ori- 
gins, they can easily be combined. For instance one easily 
checks that by carrying out a a x flip of polarization af- 
ter each spin period one cancels both linear and circular 
birefringence. 

Imperfections. In the above calculations we have 
supposed that the unknown polarization rotation B does 
not change along the fiber. The method stays valid if 
B(l) changes slowly along the fiber. To see this suppose 
that B(l) is a smooth function which can be expanded in 
Taylor series 



B(l) = B(0) + WiB(0) + O(l 2 



(16) 



We define the length over which the birefringence changes 
by 



J chge 



\diB(0)\ 



(17) 



Because B(l) is now position dependent, the evolution 
eq. must be replaced by 



17(0 - e~ iS ^ m [l-i J dire+ il ' 8 ^ ff d v B{0)ae- il ' 8( - ^ 



_ e ~iB(0)al 



I + 0{l 2 diB)) 



(18) 



This expression must be inserted into eq. ijHJ ■ The terms 
proportional to l 2 diB(0) will not cancel, so after the con- 
trol sequence, one will be left with the evolution 

U(l seq ) = I + 0(B 2 l 2 seq )+0(l 2 seq diB) 

= I + 0(u°) + O^BxBol 2 ^) + 0{uBil 2 eq /L chge ) 

where we have kept the terms of first order in cj, and we 
have supposed that diB\jB\ ~ diB /B = L~ hge . Note 
that in modern fibers with low birefringence the beat 
length Lb is comparable to the length L c h ge over which 
the birefringence changes. Hence the two terms in eq. 
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(|19fl are of comparable magnitude and there may not be 
much to gain by going to a sequence such as eq. 
which cancels the PMD to second order 

The second source of imperfections we consider is the 
fact that in general the controlled polarisation rotations 
will not be achromatic, and will therefore themselves in- 
troduce PMD. Thus, instead of carrying out the rota- 
tion o~ x , one carries out the rotation exp[icr :E 7r(l + w/3)/2] 
where (3 ~ Bi/Bq ~ In order for the method to be 
effective we must cancel the PMD due to (3. To do this 
we propose using the idea of stacking several chromatic 
waveplates to yield an essentially achromatic waveplate 
[lflj . For instance the following stack of 3 n rotations, 
with C = 1/2 and S = v/3/2, cancels all PMD to order 

e ia m (TT+uJl3)/2 e -i(C<7 x -Say)(TV+ujl3)/2 e ia m (TT+uJl3)/2 

= i(Ca x + Sa y ) + 0(lu(3) 2 . (19) 

Thus, by replacing each controlled polarisation rotation 
by three controlled polarisation rotations, one can create 
achromatic sequences. 

A third limitation is that the controlled polarisation 
rotations will differ slightly from 7r rotations in a random 
way. Thus instead of implementing o~ x one implements 
exp[i(7rcr x + r'-<7)/2] = ia x (l + if-(7 + 0(r 2 )) where r and 
f 1 are small random vector (which are simply related to 
each other). The main effect of these errors is that the 
stacks of 3 7r rotations mentioned above will no longer 
be perfectly achromatic. If we suppose that the random 
vectors r, associated with the successive controlled polar- 
isation rotation are independent, then we have to include 
an error term of the form 0(ujf3V Nr) where N is the to- 
tal number of controlled polarisation rotations and r is 
the average length of the random deviations r,-. 

As an illustration of the effects of imperfections let us 
consider the case where the residual PMD is characterised 
by length scales Lb — L c h ge ~ 10w. We correct for the 
PMD using the simple sequence of polarisation rotations 
cq. (|10|> . and in order to cancel the PMD of the polari- 
sation rotations themselves, we replace each rotation by 
a triplet of rotations as described above. After one beat 
length the Jones matrix is of order 

U(L B ) = 1 + 0(cj°) + OiuBxl^) + 0((3ujr^L B /l seq ) . 

Assuming r = 1 /30, and inserting omitted numerical con- 
stants, we find that the PMD is reduced by a factor of 5 
for the optimal value of l seq = Lb/7. 

Conclusion. We have presented a method to reduce 
PMD in optical fibers based on the use of controlled po- 
larisation rotations. The results presented here are only 
a first overview of the potentialities of the method, and 



we expect that significant improvements of the sequences 
presented can be achieved by using the same methods 
which are succesfully used to optimize pulse sequences in 
NMR. The method presented here can be combined with 
the well established method of spinning the fiber. The 
efficiency of the combined method will depend on the 
detailed properties of the residual birefringence in spun 
fibers, and on the precision with which the controlled po- 
larisation rotations can be implemented. In addition to 
the important application for long distance, high speed 
telecommunication in optical fibers, this method may 
find applications in other systems in which one wants 
to reduce unwanted birefringence. Finally, on the con- 
ceptual side, our work provides a simple system in which 
to test the ideas of bang-bang control. 
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and by the European Union projects RESQ and QAP. 
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